Magnetic reconstruction at (001) CaMnOa surface 
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The Mn-terminated (001) surface of the stable anti-ferromagnetic insulating phase of cubic 
perovskite CaMn03 is found to undergo a magnetic reconstruction consisting on a spin-flip process 
at surface: each Mn spin at the surface flips to pair with that of Mn in the subsurface layer. In 
spite of very little Mn-O charge transfer at surface, the surface behavior is driven by the e g states 
due to d xy — > d z 2 charge redistribution. These results, based on local spin density theory, give a 
double exchange like coupling that is driven by e 9 character, not additional charge, and may have 
relevance to CMR materials. 
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Despite the abundance of work on manganese-based 
perovskites in the attempt to understand the rich 
panorama of their bulk properties, 0J^] very little is 
known about their surfaces JIH or interfaces. The phys- 
ical mechanism inducing the so-called colossal magne- 
toresistance (CMR) in La 1 _ a .D, J .Mn03 (with D a diva- 
lent alkaline earth ion, and x ~ 0.3) is yet to be fully 
understood, although it seems clear that the almost half- 
metallic nature (i.e. the complete spin-polarization of 
the electrons at the Fermi level) plays a decisive role. The 
clearest evidence so far of half-metallicity is from pho- 
toemission spectra, j| whose surface sensitivity makes it 
essential to know the electronic structure of the surface it- 
self. The surface introduces the likelihood of square pyra- 
midal coordinated Mn, which is also essential to the un- 
derstanding of oxygen-deficient perovskite manganites. 
Q In addition, the interfacial behavior that is critical 
in producing low field CMR in polycrystalline material 
|jj and trilayer junctions || will involve closely related 
effects due to symmetry lowering of the Mn ion. Thus, 
studies of the surfaces of manganites are timely. 

The electronic structure of the stable phase of 
bulk CaMnC>3 has been actively investigated in re- 
cent years |]^||- CaMnOs is a G-type antiferromag- 
netic (AFM) semiconductor. The nominal ionic pic- 
ture Ca 2+ Mn 4+ C>3~ with spherical Mn d 3 configuration, 
makes the cubic (fee) phase stable jl(| over possible dis- 
tortions observed for instance in LaMnOs. In the G-type 
arrangement all nearest neighbors in the simple-cubic 
sublattice of Mn have spin-antiparallel orientation. The 
chemical picture of Mn 4+ ions is represented by com- 
pletely occupied d l\ g states. An energy gap of ~ 0.4 eV 
separate them from the empty d ej orbitals. Hybridiza- 
tion with O p states reduces the 3 [Ib nominal magne- 
tization of Mn to ~ 2.5 /is, whereas magnetic moments 
on O or Ca are zero by symmetry. 

In this paper we study the simplest manganite sur- 
face, the (001) surface of cubic CaMn03, to determine 
the surface-induced changes of structural, electronic, and 
magnetic properties. We find unexpectedly rich effects of 
surface symmetry lowering: a spin-flip occurs on the sur- 



face Mn ions that can be traced to surface states that re- 
distribute charge and spin amongst the various Mn d sub- 
orbitals and render it metallic without doping. The net 
effect is a short-range double-exchange-like phenomenon 
that relates metallicity and spin alignment, analogous to 
the CMR phases. 

Calculations were done in a local-spin-density frame- 
work; the exchange-correlation potential formula by 
Perdew and Zunger |llj was used. A plane- wave basis 
with 30 Ryd cut-off energy and Vanderbilt pseudopo- 
tentials |12| make the computation viable. To establish 
the accuracy of our methods, which have only recently 
been applied to magnetic materials, in Table || we report 
our results for the bulk CaMn03 in different magnetic 
phases. As already shown, for manganese perovskites 
the local spin density approximation successfully predicts 
the observed stable phase not only against the strongly 
unfavoured paramagnetic (PM) phase, but also in com- 
petition with closer configurations like the ferromagnetic 
(FM) and the A-type AFM (made of (001) FM layers al- 
ternating along [001] direction). Our results are in very 
good agreement with those of previous all-electron linear- 
augmented plane- wave calculations Also, our calcu- 
lated value for the equilibrium lattice constant of G-type 
AFM phase (3.735 A) is in almost perfect agreement with 
the experimental value (3.729 A). 

The stacking along [0011 consists of alternating Mn02 
and CaO units (see Fig. [l]), and the surface unit cell of 
G-type AFM CaMn0 3 is V2 x spl with respect to that of 
the bulk cubic cell (we neglect the very small structural 
distortion 0). The (001) layers are individually AFM 
and neutral, so the surface is formally non-polar. Sur- 
face formation produces two different surfaces, i.e. Mn- 
tcrminated and Ca-terminated. The presence of two in- 
equivalent surfaces in a slab will produce fictitious fields 
in the vacuum that could affect the electronic and mag- 
netic structure at the surface. We are interested on the 
Mn-terminated surface, since on it the effects on mag- 
netic properties due to the surface formation are most 
visible. Thus we use a slab containing two identical Mn- 
terminated surfaces, with mirror symmetry in the central 
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Mn layer (in total a 46-atom slab with 9 layers of atoms 
and 3 of vacuum) . 

Surface neutrality generally favours the stability of 
the ideal surface against reconstructions involving strong 
changes of symmetry and atomic density at surface. 
Therefore in this work we consider the surface with re- 
laxed but structurally unreconstructed structure, with 
different types of magnetic order. Thus we will speak of 
'reconstruction' in purely magnetic sense: on the unre- 
constructed surface the spins are oriented as in the bulk 
(Figure |l|) , while the reconstructions involve spin-flips on 
the surface layer. 

The structure of the two configurations that can be 
obtained by flipping surface spins are pictured in Fig. ^. 
In the left panel all surface spins are nipped, so vectors 
(±a, ±a) remain AFM translations but each surface spin 
is aligned with its subsurface neighbor (spin-flip AFM: 
sf-AFM). In the right panel only one (of two) surface 
spins is flipped, leaving a FM surface layer (spin-flip FM: 
sf-FM). 

Magnetic and relaxation energies and workfunctions 
for the three phases are reported in Table [fij. The AJEJ's 
reported in Table |l| are the energies gained by relaxing all 
the atoms into the slab from their ideal positions. They 
are small and reflect the small inward atomic displace- 
ments ~ 1% of the cubic lattice constant. This indicates 
a low excess stress due to the surface formation, and sug- 
gests that structural reconstructions are unlikely. The 
workfunction depends very little on the spin arrangement 
but is largest for the most stable surface. Most signifi- 
cantly, the sf-AFM surface is stable against the unrecon- 
structed one, whereas the sf-FM is the most unfavoured. 
Thus, a quite intriguing physical picture follows: at the 
surface each spin prefers to pair with the one in the sub- 
surface layer, while still keeping the AFM arrangement 
in-plane. 

It is possible to express the energy differences for dif- 
fering types of magnetic order in terms of exchange con- 
stants in a Heisenberg model 

H = — JijSi ■ Sj (1) 

<ij> 

where Sj is a unit vector in the direction of the moment 
on site j, and the sum is over distinct pairs. For the bulk 
we get first and second neighbor constants J\ = —26 
meV, J2 = —4 meV. This small value of J2 suggest that 
the nearest neighbor (nn) exchange constants contain the 
important contributions. From surface energies we get 
the nn coupling parallel and normal to the surface of 
Jll= -22 meV, J 1 - = 29 meV. While J 11 is close to the 
bulk value, has the opposite sign and is larger in 
magnitude, indicating the FM alignment of surface and 
subsurface spins is robust. 

The reversal of the surface-subsurface coupling can be 
traced to redistribution of d suborbital occupations com- 



pared to the bulk, due to the occurrence of surface states. 
The orbital-projected Mn Ad density of states (DOS) of 
the stable sf-AFM phase near the Fermi level, shown in 
Fig. [| makes evident the surface states that lie within 
the bulk band gap, which extends from -0.3 eV to 0.1 
eV relative to Ep. The surface states are of two distinct 
types, d z 2 and d xy , reflecting the strong symmetry low- 
ering of both e g — > d z 2,d x 2_ y 2 and t 2g — > d xy , {d xz ,d yz ) 
manifolds. The surface states are almost completely po- 
larized, a result of the large spin splitting A ex —2 eV 
that strongly inhibits hopping between ions of different 
spin. 

Figure |] presents the surface band structure, where for 
clarity, only the energy region of interest (roughly span- 
ning the bulk gap) is shown. The d xy and d z 2 surface 
states are easily identifiable. The d xy states at the sur- 
face are shifted upward by 1 eV and overlap with d z 2 
states that in the bulk hybridize strongly with the O pa 
orbitals and form low-lying bonding and high-lying an- 
tibonding bands (most of the d weight is in the latter). 
The two d z 2 states (one from each surface of the slab) 
are split by 0.2 eV as a result of the interaction between 
Mn at opposite sides of the slab. (For a thicker slab they 
would converge to a single, doubly degenerate band, av- 
eraging the calculated d z 2 states). The d xy band has a 
bandwidth of 1.4 eV and a dispersion that follows 

e 3 ^ — — 2t[cos{k x + k y )a + cos(k x — k y )a] (2) 

using the conventional perovskite coordinates. This dis- 
persion arises from hopping between second nn, which 
are the nearest neighbors of like spin. The effective hop- 
ping amplitude is t = 0.17 eV. The d z 2 band is very 
narrow (0.2 eV) and its dispersion is not easily repre- 
sented by tight binding form, reflecting small competing 
hopping processes along the surface (and perhaps sub- 
surface) that are not easily identified. Coupling of the 
d z 2 state perpendicular to the surface is large, however, 
as reflected in the penetration of the state onto the fifth 
atomic layer (third Mn layer). 

The net effect on the Mn ion of surface formation is 
an intra-atomic shift of charge from d xy to d z 2 orbital, 
without appreciable change of charge or moment. In the 
solid the magnetic moment comes almost entirely from 
t2 g states, whereas on the surface d z 2 (surface) states con- 
tribute about 30% of the moment. Only the subsurface 
Mn in the sf-AFM phase experiences a net gain (0.07 hb) 
due to the partially occupied d z 2 state not compensated 
by a depletion of d xy states. In bulk, magnetic moments 
in the G-type bulk are only allowed on Mn, by symmetry. 
With the surface formation, O in-plane with Ca acquire 
a magnetic moment as well. This is larger in the sf-AFM 
(0.11 1-J.b) than in the unreconstructed phase (0.06 /^b), 
because it is enhanced by the parallel alignment of two 
neighboring Mn spins. 

Most of the characteristics inferred by DOS and band 
structure analysis can be better visualized by means of 
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isosurface plots (Fig. ||) of charge density and magneti- 
zation of the stable sf-AFM phase. The quantities shown 
are due only to states within the bulk gap (see Figure ||), 
thus represent charge and magnetization of the surface 
states. The charge density clearly shows both d xy and 
d z 2 characters of the charge on Mn, as well as a p^-type 
contribution from O. 

The physical mechanism driving the changes in ex- 
change interaction parameters at the surface is related to 
that described by Solovyev et al. |l3| investigating how 
the Jahn- Teller distortions (JTD) affect the magnetic or- 
dering of LaMnC>3. The basic driving force in the FM- 
to-AFM transition of LaMnC>3 vs. JTD is the decreasing 
of d z 2 occupancy occurring with JTD. The d z 2 — d z 2 in- 
teraction is indeed a dominant positive (i.e. FM) contri- 
bution to J . Also positive are the d x 2_ y 2 — d z 2 and the 
much weaker d x 2_ y 2 — d x 2_ y 2 interactions, whereas t2 9 
orbitals interact by superexchange, and their contribu- 
tion is AFM 1 1 3 1 . When the d z 2 orbitals are sufficiently 
occupied to make J 1 - larger than the next nearest neigh- 
bor interactions (favouring AFM) the order becomes FM 
along the z axis. 

In the present case the picture follows analogously: 
surface formation, not JTD, results in dehybridization 
and partial filling of the d z 2 states on surface (and sub- 
surface) Mn, and a partial depletion of d xy orbitals. As 
a consequence J 1 - changes sign and the magnetic order- 
ing along z reverses. A verification of this mechanism 
is given by comparison of the band structures (or DOS) 
of the two competing phases (Figure ^): in the sf-AFM 
phase there is more d z 2 occupation and less d xy depletion 
than in the sf-FM phase. Also, the slight occupation of 
d z 2 states on the subsurface atom is larger for the recon- 
structed sf-AFM phase, but not sufficient to propagate 
the spin alignment further into the bulk. The consid- 
erable difference with respect to the LaMnC>3 JTD is 
that its distortion is extended, whereas at the surface 
the ordering is a local effect limited to the first two lay- 
ers. This local spin-flip process is likely to be relevant 
to more general situations in manganites, such as at sur- 
faces and interfaces of doped systems where it may affect 
spin transport, and at Mn sites neighboring O vacancies, 
as in CaMnOs-z. [| 

To summarize, we have described a spin-flip process 
at the Mn-terminated (001) surface of CaMn03 that is 
driven by symmetry lowering due to surface formation 
which causes the partial occupation of the e g d z 2 surface 
states. This partially occupied narrow d z 2 band may dis- 
play correlated electron behavior. This d z 2 occupation 
reverses the magnetic alignment (from AFM to FM) at 
the surface in the direction orthogonal to the surface but 
conserves the AFM symmetry along the surface. The 
surface states are almost completely polarized, but AFM 
symmetry requires that both spin states occur in equal 
number, so this result may be difficult to verify experi- 
mentally. 
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PM 


FM 


A AFM 


G AFM 


E(eV) 


+0.56 


+0.16 


+0.08 





M Qu s ) 





2.56 


2.48 


2.36 



TABLE I. Energies (per formula unit, referred to that of 
the stable G-type AFM phase) and magnetic moments on Mn 
atoms for the bulk CaMnOs in different magnetic phases. 





UNREC 


SPIN-FLIP FM 


SPIN-FLIP AFM 


E(eV) 





+0.12 


-0.12 


AE(eV) 


0.09 




0.05 


W(eV) 


5.65 


5.64 


5.73 



TABLE II. Energies, per surface V2x^2 cell, of the (001) 
CaMnOs surface in different magnetic phases (see the text); 
AE'(eV) are the corresponding relaxation energies, and W the 
workfunctions. 
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FIG. 1. Structure of the Mn-terminated (001) surface of 
G-type AFM CaMnOs in the unreconstructed magnetic con- 
figuration, i.e. the spin orientation at surface is equal to that 
in the bulk. 
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FIG. 2. Two possible magnetic reconstructions of 
Mn-terminated (001) surface. Left picture: spins on both Mn 
at surface are flipped, the surface is still AFM. Right picture: 
spin on one of two Mn is flipped: FM surface. 
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FIG. 3. Orbital-projected DOS of Mn d states for the (001) 
spin-flip AFM surface. Upper, middle and lower panels re- 
fer to surface, sub-surface and central Mn atoms in the slab, 
respectively. Dashed, solid, and thick solid lines refer to 
d xz +d yz , d z 2, and d xy orbitals, respectively. Up spin DOS 
is plotted upward; there is another atom in each layer whose 
DOS is flipped from those shown. 
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FIG. 4. Band structure in a small energy window around 
the Fermi energy. Left panel refers to the unreconstructed 
surface, right panel to the spin-flip AFM. Dashed lines 
are boundaries of the bulk energy gap. High symme- 
try points in the IBZ are F=(0,0,0), X=(27r/a')(l/2,0,0), 
M=(27r/a')(l/2, 1/2,0), in the y/2 x ^2 surface cell. 




FIG. 5. Isosurface plots of the surface states for the 
spin-flip AFM phase. Left panel: top view of the elec- 
tron charge isosurface (of 0.01 electron/Bohr -3 magnitude) 
at (100) surface. Right panel: tri-dimensional view of the 
magnetization. By inversion symmetry along z only half slab 
(five layers) is shown; on top there is the vacuum, on bottom 
the bulk. Light and dark isosurfaces are of same magnitude 
(0.005 /is/Bohr -3 ) and opposite sign. 



